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Наближений анал!з процесу кристал!1зацИ металу 


Рассматривается задача тепловой обработки металла. Строится приближенное решение задачи. Получены 
оценки для температуры и теплового потока. В ходе численных экспериментов определены параметры 
тепловых потоков при отделении от стенок кристаллизатора. 
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Тье ргоет оЁ Феппа! ргосеззтю оЁ теа| 15 сопз14еге4. Тве арргохипае зойоп 15 сопзгасе4. 
ЕзитаНопз Юг а {етрегайиге ап {Пегта]| эгеат аге 201. ш фе питенса| ехрегипеп&5, е рагатейегз оРе 
Беаё Но\/ тю зерагае Нот йе \уаП$ оРе то14. 


Кеу могд$: шооб, сгуаШиег, те], птитииайоп ое РпсНопа|, е!есёгоз1ас тете то. 


Розглядаеться задача теплово! обробки металу. Будуеться наближений розв’язок задач1. Доведен! 
ошнки на температуру 1 тепловий потк. У ход! чисельних експериментв визначено параметри 
теплових поток при в1ддленн! вд стнок кристал!затора. 


Ключов: слова: злиток, кристал1затор, метал, мим1защя функцюналу, 
електрошлаковий переплав. 


Распространение тепла в различных средах оказывает большое влияние на характер 
протекания многих важных для практики процессов. Среди задач, связанных с распро- 
странением тепла, выделяется класс задач, в которых исследуемое вещество переходит 
из одной фазы в другую с выделением или поглощением тепла. 

Целью данной работы является моделирование процесса кристаллизации метал- 
ла, изучение процесса завершения получения слитка в кристаллизаторе путем его вы- 
тягивания. 

Рассматривается задача управления информационными процессами при авто- 
матизации технологий тепловой обработки металла, на основе математического мо- 
делирования, анализа статистических данных и теплофизических экспериментальных 
измерений. В качестве источника информации исследуется математическая модель, 
основанная на пространственной задаче Стефана, с учетом конвективного движения 
и примесей в жидкой фазе. 
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Постановка задачи 


В настоящее время электронно-лучевой переплав (ЭЛПИ) является одним из 
наиболее эффективных способов повышения служебных характеристик металлов и 
сплавов. Он применяется в исследовательской практике и промышленности для по- 
лучения особо чистых металлов и сплавов с высокими физико-химическими и меха- 
ническими свойствами. 

Накопленный к настоящему времени экспериментальный материал по опытно- 
промышленному применению электронно-лучевого переплава не только подтверждает 
его преимущество по сравнению с другими методами вакуумной металлургии, но и 
выявляет ряд преимуществ. ЭЛП является наиболее приемлемым и экономически 
оправданным процессом получения особо чистых ниобия, тантала, циркония, вана- 
дия, меди, никеля, железа и других металлов и сплавов на их основе. Исключение 
составляют металлы с высокой упругостью паров в точке плавления, например хром 
и сплавы на его основе. 

Весьма эффективного применения электронно-лучевого переплава для повыше- 
ния качества специальных сталей и сложнолегированных сплавов на основе никеля и 
железа. При переплаве указанных материалов одновременно с резким уменьшением 
содержания водорода, азота и кислорода происходит практически полное удаление 
многих легкоплавких примесей (свинца, цинка, висмута, олова и др.), снижается 
содержание неметаллических включений. Наблюдается также диспергирование и более 
равномерное распределение по сечению слитка оставшихся в металле неметаллических 
включений. 

Значительный прогресс в области электронно-лучевой технологии достигнут в 
результате перехода от прямого переплава расходуемой заготовки в кристаллизатор 
к переплаву с использованием промежуточной емкости (ЭЛПЕ) и способу формиро- 
вания слитка в горизонтальном кристаллизаторе (ЭЛПГ) [1]. 

Пусть Б=(-1<х<|1у<0) — полуполоса, заполненная твердым металлом. 


Обозначим через и(х, у) температуру этого металла. В работе [2] построена функция 
и(х, у) как решение соответствующей краевой задачи. 
Отождествим теперь температуру и(х, у) с температурой твердого слитка, на- 


ходящегося в кристаллизаторе при электрошлаковом переплаве. Для вытягивания 
слитка из кристаллизатора поверхность слитка предварительно обогревается тремя 


электронными лучами Й”,И’, и И’, ‚ причем мощность И’, одного из них равномерно 
распределена в центральной зоне 1 <х>Ьу= 0} ‚ а два других сконцентрированы 
по краям х = +1 [2]. Независимо от того, в каком отношении находится температура 


поверхности слитка с критической температурой Т”, при которой поверхность слитка 
отделяется от стенок кристаллизатора, теплообмен слитка с кристаллизатором осущест- 
вляется по формуле 


и, аи = 0,х = +1-® <у<0. 


Для получения температуры слитка достаточно положить у(х) = (И, И’,,И’,). 
Далее, введем в рассмотрение функционал 


И) = [(м.у)-Т” 4. 
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Рассматривается задача. Требуется определить поток у(х) из допустимого 
множества И, доставляющий наименьшее значение функционалу Г(у). Минимизирую- 
щая последовательность у, строится по формуле у. =У, +=, (у, -у,), параметр 
=, выбирается из условия шш /(у, +=, (у, -,,)), О<&, <1 [3]. В качестве области 
определения функции И берется множество кусочно-постоянных ступенчатых функций: 


УЕУ,,Х, хх, ,У, = СОП5Е, К = 0/1,2,,...„т. 


При этом формула (5) примет вид 


С с05 И„хе“”” эшЯ,х, -зШ Ах, 
и(х, У) —= У. 2 У. У, > 
и с0$° Я, ко Я, 
И» С р 0 12 ) 


а Г(у) = Ку, уу, ,...Уи). 

При численной реализации задачи необходимо учесть ограничение [4] 
2500 <у(х) <5000, здесь у(х) — мощность потока в единицах МВт/м”, а также 
@ = 2,66 (число Пекле), @ = 3,05 (число Нуссельта). 


Способы решения задачи 


Построение минимума функционала /(у), 
0 


| 2соз хе”  зшд 
где Ку) = [(иру)-Т”)? у, иобк,у) = Ав = бурь, 


с05°Ж. А 
и ибо“ 
М 
@ |9 о) _ 
др = = + тт +1, п=1.2,3...1, — положительные корни уравнения А = @,с14. 


Обозначим через и(х,у) температуру металла. Отождествим теперь темпера- 
туру и(х, у) с температурой твердого слитка, находящегося в кристаллизаторе при 
электрошлаковом переплаве. 


Необходимо найти шш/и у,. Минимум функционала / находим из условия 
О1 й 
——=2 [а у»э-тл, у =0. 
ду 


Непосредственные вычисления показывают, что 


ет НТ”, ет 

с0$ с0$ 

2 (1+ 9 р} г. Дол (1+ @ 2 >) 
( ( 

Второе приближение. Требуется найти минимум функционала 


1) = [м 4, у)-Т")? 4. 


Ку) = 
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если 
2со5 хе“  зшА 2со5 д хе“  зшА 
и, (х,у)= — ту+ —- ти = [% (х, у) + Ло, Ур. 
с05 ЖЖ. № 05 Я. 4% 
№ 60° р} ) А+ в, р} ) 


0 1 
Необходимо найти у, и [(у,). Найдем теперь минимум функционала /, в 


случае когда и (х, у) = (Д(х, у) + Л (х, у), 
со; йхе“" ша 
где Л(х,у)=2 — 
со 
и (+9, —— ) 


1 


Поступая, аналогично тому, как это было сделано в случае нулевого прибли- 
жения, получим 


5124 (1-е^") - 5124, 1-е") 
с05? Я с05? 7 
Лоо (1+ @ 2 2 ) Л (1+ © р} ) 


0 1 


ОР 


т? %(1-е*“"”) ы эт? 24, (1-е?!) 


А= 
2 2 
с05° Я с057 Я 
Лод! А+ @, 2 °.)? 2щА 1+9, 5 —): 
где 0 1 
ыы т 2 зт 24, (1-е) 1-е") 
с05° А) соз? д 


Ло (1+ @, ой (1+ ©, р 


о 1 
Далее, имеет место следующая формула: 
ПР ее" ш24, (1-е^" 
аа $12 (1-е — +2 $ш24, (1-е — р 
322 05° М. 2 372 с0$° 7, 2 
2 + @, 2 ) 2 А +9, 2 ) 


о 1 


5124, зт 2 (1-е 1-е") 
с05* А со$° 7 
Лоо (1+ у 2 ° Уи + @% р} -) 


( 1 


+ 2, 


у ш2%, (1-е^" 124, (1-е“й Е 
т" $12, (1-е > эм ‚1-е _ НТ? 
с05 с05 
Я (1+ @ р} °) Яд (1+ @ 2 ) 


( 1 
Минимизация /[(у) при ступенчатой функции у(х). 


1 


ет "66034,54, 


п\0 
р 
7 


здесь у(х) =у, при 0 <х<х,, у(х) =ъ, при х, <х<1, М, < У(х) < М.. 


К) й ИУ 
Имеем и(х, у) = 2. А - 


= И, (1 т 5 
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В дальнейшем положим у, = №, и К (х, у) =2 


( 
0 
. 7 
Требуется найти шш/Ги у, ‚ где Г(и,М№,) = | (1, У(аи, +ЬМ№,)-Т’)* 4», т.е. = =0. 
Н и 


Имеем 
—1 


‚ с084(-е^”) со5? (1-е?) т со5? 41-е?” 


и =|27Т г) —25М№, 7 1 ‚ 
с05 с05 с0$ 
шо (1+ ©. 2 °) (1+. 2 °) (1+ 6% 2 ° 
Жо Ао о 
а 5ш Йох, и $ш Я, — 51 Фух, 
Ло № 
0 
Ци, М№,) = |, Уау, +Ьм,)-Т`)? а = 
НЕЁ. 
2 
26087 1-е?" : %-е^# „2 
Е о 
2 с05° № Ио 2 с05° 
№ 1+ @ 5 ) Шо (1+ @% 5 ) 
№ Л 
Далее, имеем 
й У 1 А хе“ 1 
м = [ив ооз ль ваЯ +2 АН [и()соз ва = 
с0$° о А 
(+ @ 5 ) и + 9 =) 
ие И 
Я. оу Я. МУ 
Не ПЕ, 
пана с0$ №. а с0$ ) 
1 0 
Е У Яух, = Ш Я, — за дух, а эт ИХ, = шт, -зш Ах, 
Ло Ло И И 
с0$ Ди хе“ с0$ хе“ 
№ (ху) =2 о 2 ‚ (ху) =2 2 
с05° № с0$° #, 
(1+ @, р} ) д + 6, р} ) 


( 1 


Таблица 1 — Значения тепловых потоков 

(а № м т Н Г (00) Го) | йо) 
0,4 11,2646 | 3,9352 | 0,9 —3,0 | 226,77 | 148,999 | 59,9 
0,4 11,2646 | 3,9352 | 0,9 —4,0 | 237,580 | 58,868 75.6 
0,4 | 1,2646 | 3,9352 | 0,9 —6.0 | 241,625 | 166,516 | 95,1 
0,4 11,2646 | 3,9352 | 0,95 | -3,0 | 244,449 | 147191 | 66,7 
0,4 | 1,2646 | 3,9352 | 0,95 | 4,0 | 264,854 | 156,838 | 84,2 
0,4 | 1,2646 | 3,9352 | 0,95 | -6,0 | 238,346 | 164,304 | 100,5 
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Были проделаны численные эксперименты при определенных параметрах [5]. 
Полученные результаты были представлены в табл. 1. Анализ численных результатов 
показывает убывание значения функционала на первом и втором приближениях, что 
соответствует смыслу теплофизических процессов происходящих при электрошлаковом 
переплаве (ЭППУ. 


Список литературы 


1. Патон Б.Е. Избранные труды / Патон Б.Е. — Киев : Институт электросварки им. Е.О. Патона 
НАН Украины, 2008. — 893 с. 

2. Шевченко А.И. Моделирование одного класса сложных систем с нечетким управлением / А.И. Шевченко, 
А.С. Миненко, И.А. Сыпко // Доп. НАН Укра1ни. - 2013. — № 8. - С. 52-54. 

3. Сыпко А.И. Приближенное моделирование процесса кристаллизации при наличии конвекции / 
А.И. Сыпко // Искусственный интеллект — 2013. — № 2. - С. 80-85. 

4. Сыпко А.И. Численное моделирование процесса кристаллизации металла / А.И. Сыпко // Искусственный 
интеллект — 2014. — № 1.— С. 153-159. 

5. Будак Б.М. Сборник задач по математической физике / Б.М. Будак, А.А. Самарский, А.Н. Тихонов — 
Москва : Наука, 1980. — 686 с. 


ВеЁегепсе$ 


1. Раюп В.Е. Зеецеа угогКз / Раюп В.Е. - Клеу: Иесёлс \Уе!ат® ш5ийие патеа Е.О. Раюп оЁР фе МАЗ 
оЁОКгаше, 2008. — 893р. 

2. ЭЗКеусвепКо А.1., МшепКо А.5., Г.А. Зурко Мо4еПп» о фе опе с!аз$ сотр! ех зузетз мВ Ё172ту 
сопйг`о] // Керог$ о МаНопа| Асадету о{ Зс1епсез о ОКгаше —2013. — № 8. -Р. 52-54. 

3. ЗурКо 1.А. Арргохипае то4деПп» оР фе сгузбаШтайноп ргосез$ ш Пе ргезепсе оЁ сопуесйоп // Агйс1а1 
шеШоепсе —2013. — № 2. - Р. 80-85. 

4. ЗурКо 1.А. Матепса! поете оР@е сгузбаШтаноп ргосезз оЁтеа! // Агийса! пеШоепсе. 2014. № 1.Р. 153-159. 

5. Виаак В.М., ЗатагзЮу А.А., ТБопоу А.М. СоЙесноп оЁргоепаз оп тайетайса! рпуз!1сз. — Мозсо\: Заепсе, 
1980. — 686 р. 


ВЕ5ОМЕ 
Т.А. ЗурКо 


М№Митетса! Модейие ог те СгуяюШ2аноп Ргосез; ор Ме 

ВасКогоипа. Атопо фе ргоептз аззослае4 ул Фе Ч151ийоп оЁ Беаф, Фете 15 а 
с1а5$ оЁ ргоетз ш УеЬ фе шуезиае4 заб$апсе раззез Нот опе рВазе 0 апоШег, 
У Ше гееазе ог абзогрНоп оЁВеай. 

Мажена[5 ап шео9$. ТЬе рогрозе оЁ {1$ уогК 1$ ю забапйае Фе тае- 
тайса| то4е| оГ Фе ргосез$ оГ сгубаШтаНоп оР теа|. ТВе ргоет оЁ сопёго| оЁ 
{есппо|ос1са| ргосез$ оР те{а15 Фегта| ргосеззто 1$ сопз14еге4. Аз шЮоптайоп гезопгсе 
Ще Шгее дитеп$1оп сопуесйоп ЗеЁРап ргоМет ш Паи1А рБазе 15 шуезиоае4. 

Везий. ТЫ$ рарег ежепа ю ите-дерепдепй сазе зоте гезаЙ оМашед Бу Фе аифог Юг 
уеаду-5ае ЗеРп ргоМет уз сопуесйоп. То деегтше Фе оршта! Фета! соп4 010$ оЁ 
Ююгтайоп оЁ 11201 сасшаНопз \уеге саглей ой ш Фе Нате\уогК оЁ фе тафетанса| поде! оЁ 
Фегиа| ргосеззез ш а суппанса! 120%, адар4 ог Фе сазе оЁ БоПо\и шеоб. ш фе тоде| оЁ 
Подаа те 15 роге шю Фе той 9 рогНопз, ап4 а \едее оЁ И 1$ риЦе4 ой репо саПу. 

Сопсм$1 оп. А тафетайса| под ус аПо\з ю амепише е Фегта| поде ш уфись 
[1$ дебпед. Зиича 4 ргосез$ оЁ сгузаШтаноп оЁ теа|, раззште ш зрес1а! теаПагоу, патейу 
эфиед Пе ргосезз оРсотреноп оРШе тгесе!ре оРап тео ш Фе тошА Бу змеюНЫто. 

Тве рго ет оЁ Фегта| ргосеззште, оР тейа| 15 сопз14еге4. ТВе арргохитае зошйоп 
1$ сопзнаце4. ЕзИтаНопз$ г а {етрегайге ап Фегтла| згеат аге 201. 

ТБе г@еуапсе оГ Фе ргезеще4 уго!К 1$ дие Бо 10 Ше ргасИса! детап4$ ог №77у 
соп@го| оЁ Пе ргосез$$ оЁ сгубаШтаНоп ог Фе оБ]ес{ ул сотр]ех оеотецу. 
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